ON THE GAPS BETWEEN ZEROS OF TRIGONOMETRIC 

POLYNOMIALS 
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Abstract. We show that for every finite set ^ S C Z d with the property 
—S = S, every real trigonometric polynomial / on the d dimensional torus 
T d = R d /Z d with spectrum in S has a zero in every closed ball of diameter 
D (S), where 

and investigate tightness in some special cases. 



1. Introduction and presentation of the main results 

We are interested in the following question: for real trigonometric polynomials 
with a given spectrum, what is the maximal possible distance between two con- 
secutive real zeros? Of course, if we can use enough frequencies for building our 
trigonometric polynomial, the size of the largest gap can be as close to 2ir (or to 
1, if one uses the transformation used throughout this paper) as we want. The 
question gets more interesting when the number of frequencies is relatively small. 

The first theorem is a very general result, partially answering the above question. 
We state it for multivariate trigonometric polynomials on the d-dimensional torus 
T d = R d /Z d . 

Theorem 1. Let S C 2, d be a finite set s.t. —S = S. Let 

f( x ) = ^2 C ( A ) exp 27r i(x, A) 

be a real value trigonometric polynomial on T d . Then f has a zero in any closed 
ball of diameter D(S) where 



D(s,: =£w 



Note that the requirement — S — S is necessary because we are interested in real 
trigonometric polynomials. 

How sharp is this theorem? To answer that we have to analyze specific cases. 
From now on we concentrate on the one dimensional situation. For a finite set 
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O^ScZ define 

M (S) := sup sup |/| , 

spec /CS/C [0,1] 

where spec / is the spectrum of / i.e. the support of the the Fourier transform /. 
The second supremum is over all intervals I such that / is never zero on I. Here 
and also later on |/| stands for the length of I. With this notation theorem 1 is 
shortened to the formula 

M(S) < D{S) . 

The easiest case to investigate is when the spectrum is a (discrete) interval, or 
rather, two intervals situated symmetrically around 0. It turns out that in this case 
it is possible to get an explicit formula for the supremum (maximum does not exist) 
of the gap size: 

Theorem 2. For S = [-N - K, -N] U [N, N + if] C Z with N, K e N, 

- m< • 

As the proof will show, we can also demonstrate an explicit formula for the 
trigonometric polynomial corresponding to the supremum (being a polynomial with 
"touching zeros" — e.g. non-negative on an interval of length (K + 1) / (2N + if)). 
We see that the general result applied for this specific case is sharp only asymptot- 
ically when if <C N. This leaves many questions. As examples, we formulate some 
of them: 

Question 1. Theorem 2 can be generalized to linear progressions with stepping b 
smaller than twice the starting element, i.e. when S — ± {N + rib : n £ 0, 1, 2, . . . K} 
for some K, b £ N, b < 2N (see Theorem 3 on pageQJ. What happens when b > 2N? 
Can an explicit formula for M (S) be found in this case too? 

Question 2. What is the maximum when the spectrum points are perfect squares, 

i.e. S= {±nYj: ? 

Theorem 1 states that in this case the maximum tends to zero as N — » oo even 
if if — > oo, a result which by itself is intriguing. However, is the asymptotic truly 
0(1/N) or is it even lower? 

Question 3. What can we say when S is a random set? A typical model is taking 
every integer in 1, . . . , N with probability r and independently and then symmetriz- 
ing. For asymptotic results we should assume t and N are related somehow, say 
t = N~ s for some < S < 1. 

Question 4. What about nets? A net of order n is a set with 2™ +1 elements, 
defined by parameters {ao,...,a n } as follows: S :— ± {ao + J27=i {A a i}}- I n 
many questions nets behave like generalized arithmetic sequences. Is it true in our 
case as well? 

We should remark that the specific case of N — 1 in theorem together with 
the extremal polynomial, has been known for a while. The question was asked in 
T65J and answered in [B84J, which actually showed that this polynomial maximizes 
the measure of the positive part. The same polynomial reappears attributed to 



Arnold in T97] and again in |GS00| (which also handles the question of positive 
coefficients). For our needs, though, the case N = 1 is not enough to make a 
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meaningful comparison between M and D since we are even more interested in 
the case when they are both small, and the proofs in these papers do not seem 
to generalize nicely. See also [ViXi for inexact multidimensional results. Finally, 
BH84, MDS99J deal with the related question of completely positive polynomials 
with some bounds on the zeroth coefficient. 

2. Proof of Theorem 1 and Theorem 2 

Proof of theorem 1. We will use induction on the cardinality of S. Due to the 
symmetry of S, we have #5 = 2k, k = 0, 1, 2, . . . The case k = is trivial as then 
/ = 0. 

In general, given k > 0, assume we have proven the statement for all sets S with 
cardinality smaller than 2k. Take a set S with cardinality 2k and a trigonometric 
polynomial / (x) = ^ AeS c (A) exp 27ri (a;, A). Arguing by contradiction, assume 
that / does not have a zero, say is strictly positive, on a d-dimensional ball of 
radius R = \D(S) centered at y = (j/i, . . . ,yd), which we denote as usual by 
B(y,R). Pick a frequency v — {y\, . . . , v£) e S, and translate / by /i :— „,, v '-rm . 

Since / is strictly positive on the ball B (y, R), the translate is strictly positive on 
the translated ball B(y — fi, R). Adding the translate to / (x) we get a new function 

}{x) :=f{x)+f{x + l i) 

which is strictly positive on the closed ball 

B{y-\ii,R-\\\n\\2) CB(y,R)r\B(y-»,R) (1) 

whose diameter is 2R — M . 

The Fourier expansion of the translated function is 

f (x + fi) = c (A) exp 2-ki (x + fj,, A) 

= ^""^ c(A) exp 2m(fj,, A) exp 27ri(x, A) 
AeS 

and thus 

/ (x) = c(A)(l + exp 2m(iA, A)) exp 2ni (x, A) . 

AGS 

Since for A = ±v we get exp2ni{ii, A) = e ±m = — 1, the spectrum Sj := spec(/) of/ 
is contained in S but definitely does not contain ±z/. Thus D(S) — 2 \\l\\ 2 — D(Sj), 
and so we find that we have a new set Sf of cardinality smaller than 2k and 
a trigonometric polynomial / with spectrum Sj which is strictly positive on the 
closed ball B(y — /i,R(S^)) of radius R{Sj) — ^D(Sj) centered at y — [i. This 
contradicts our inductive hypothesis. □ 

Remark. A similar proof can be applied to shapes different from balls. For example, 
for cubes we get the following version of theorem 1: If spec / = S then / has a zero 
in every closed cube of side length L(S) where 

i(s): ^«t ■ 
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The only proof element that needs modification is (QJ , as for cubes we have (denot- 
ing by C(y, L) the cube of side length L based at y), 

C(y, L - UmIIoo) C C(y, L) n C(y - /x, L) 

and HmHoo = IMU/IMIi- 



Proof of Theorem 2. Theorem 2 is a specific case of the following 

Theorem 3. For S = {-N - Kb, . . . , -N -b, —N} U {N, N+ b, . . . , N + Kb} with 
N, K, b £N andb < 2N we have 

M(S)= K + 1 



bK + 2N ' 
Proof. First we show 

K + 1 

M(S)< (2) 

K ' ~ bK + 2N W 

by proving that for every trigonometric polynomial / (t) = J2xeS c M e 2wlAt , if / > 

on an interval / then |/| < ■ We may assume w.l.o.g. that / = [0, a]. One 

can write / (x) = ReF(i) with F (t) := e 2 ™ Nt Q (e 2mt ) , where Q is a polynomial 

of degree Kb, that is 

Kb 

Q(z) = cIJ (*-&0 • (3) 

3=1 

The condition / > on / is equivalent to aigF 6 (— 7r/2,7r/2) on /. Since e 27rlAr * 
does a rotation of 2irNa in the positive direction as £ goes from to a, / > 
on I requires that Q (e 2 '"' 1 ') does a rotation of at least 2nNa — ir in the negative 
direction on /. In other words, if argQ is any continuous version of arg on the arc 
{e 27 ™* : t 6 /} (arg exists because f\j > gives Q ^ on the arc, so for example 
take afgQ (e 27 ™*) := Im J* ^) then a necessary condition is 

hid Q := argQ (e 27rm ) - arg Q (1) < tt - 2nNa . (4) 

Clearly the index ind Q of Q does not depend on the choice of arg (all versions of 
arg differ by a constant 2irk). As we want an upper bound for a = \I\, we want to 
maximize — ind Q. Clearly 

Kb 

ind Q = ind {z — £j) 

3=1 

Now it is obvious that if £ is inside the unit disc, then ind {z — £) is positive (draw 
a picture!) — this is not what we want. If £ is outside the disc, then ind (z — £) 
is harmonic (as a function of £) and therefore has no maxima in any open subset 
outside the disc, thus we need to investigate only the limiting behavior as £ tends 
to infinity and to the unit circle. 1 Both can be calculated explicitly. When £ — > oo 
clearly ind (z — £) — > 0. On the other hand, a simple calculation shows that 

lim md(z-Z) = h-' K ifse (°' a ) ( 5 ) 
e^e=™ V W ira if s ^ [0,a V ' 



^The fact that the supremum occurs when ^ converges to the unit circle from outside can also 
be deduced by elementary geometric arguments. 
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where the limit is taken from the outside of the unit circle. A geometric note may be 
due here: the limit lim ? ^ e 2^i 3 (arg ( e 27rM — £j — arg (1 — £)) , viewed as a function 
of a, for s € I is not continuous — it has a — ir jump discontinuity at s, this gives 
the — 7r factor in ©. 

A similar calculation shows that if £ — > 1 or £ — > e 27rw then the limit of ind (z — £) 
depends on the angle of approach, but in any case is bounded from below by 7ra — 7r. 
These three cases (£ inside the disc, outside, and approaching the boundary) totally 
give 

ind(z-£) >7r(o-l) V( £ C \ {e 2ra : s £ /} . (6) 

The case b = 1 is now immediate from and ©■ The case 6 > 1 is not much 
harder. Let t = e 27ri / b . If £ is a zero of Q then so are t£,t 2 £, . . . ,t 6_1 £ so we 
can rearrange £j such that £j+n = T £,j- In particular this gives the following 
representation using the first K zeros: 



3=1 



Now the analog of © is 



ind - C b ) > 7T (6a - 1) (7) 

which holds for a < 1/6 — we will justify this assumption later. The reasoning is 
similar: if £ is inside the unit disc so are all the conjugates and the index is positive. 
For £ outside the disc the index is harmonic in £ so we need to investigate only £ 
converging to the circle, but as a < 1/6 only one of the conjugates of £ may converge 
to the arc {e 2nis : s € (0,a)}. Summing 6 copies of © gives (0. Summing the K 
copies of O gives (with J3J) 

JsTtt (6a - 1) < ind Q < tt - 2nNa (8) 

which, again, holds for a < 1/6. As according to our conditions 2N/b > 1, for 
a = r — e with e > sufficiently small (jHJ cannot hold, thus no such / can have 
an interval of length \ — e with no zeros. This justifies our a priori assumption 
a < 1/6, and we can use JHJ to get the largest a allowed, proving (2J. 
We turn to the proof of 

M ^6#^ ■ ^ 

The first step will be to construct an extremal function / which will only satisfy 
/ > (i.e. will have zeros in the interval I — [0, a] where a = gjq^jy ) ■ The analysis 
above suggests that we try / = ReF, F (t) = e 2mNt Q (e 2mt ) , for Q satisfying ©, 
with K zeros on the arc {e 2ms : se 7} (the other zeros are the conjugates, which 
are of course also an the unit circle but not on the arc). It turns out that it is 
enough to space the zeros evenly, namely define 

K 

Q{z) :=c]l(z b -^ b ) 

3=1 

where r := e 2 ™ 1 with 77 := bK + 2 N » an< ^ c ^ s some constant such that F (0) = 
Q (1) = — i. Now argF can be calculated explicitly, where this time arg is in the 
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usual sense, i.e. a discontinuous function between — tt and it. Elementary geometry 
shows that 

Summing this for all the roots we get that on each segment (ki], (k + 1) 77) arg_F is 
a linear function with derivative n (bK + 27V). Thus, on (0, 77) it is increasing from 
— 7r/2 to 7r/2, at 77 it has a jump of — 7r, then on (77, 2rj) it is again linear with the 
same slope, so it again rises from — 7r/2 to 7r/2, then has another jump of — 7r and 
so on. This argument works till (K + 1) r\ = a. Thus we get arg F e [—ir/2, tt/2] on 
/ and therefore / > on / having its only zeros at the points {0, 77, 2rj, . . . , Krj, a}. 

Taking / (t) + f (t + e) (whose Fourier transform is also supported on S) for e 
small enough one gets a strictly positive function on any interval interior to /. This 
proves © and the theorem. □ 



3. Corollaries and remarks 
Let's compare the results of theorems 1 and 2. A simple calculation shows that 

N+K i 1/ 1 1 

°W = IT. ST = T 



n=N 



2n 2N + K 



N 3 



It is interesting to note that the true value, i.e. M(S) = 2 I n+k ls ^ ne un iQ ue 
approximation of D(S) by a rational function of order 1 with this error term. 
The following fact is a corollary of Theorem 2: 

Corollary. For every interval I C [0, 1] with \I\ < 1 there exists some constant a 
such that for all N eN there exists a real trigonometric polynomial f with 



spec / c [-aN, -N] U [N, aN] 



such that f > on I . 



Proof. Let 5:=1-\I\. Let a := 1 + 4/5. Apply Theorem 2 for N, K := [(a - 1)JVJ 
([■J denoting as usual the integer value) and e := 5/2. Then as in this case 



f K+l 
\ 2N + K 



\I\ > 



l + S{N-±) 
2N + K 



>0, 



there exists a real trigonometric polynomial f\ with spec(/i) C [—a A, —TV] U 



[A, aA 
terval 



s.t. /1 > on an interval of length bigger than |/| (namely, on the in- 



0, 



K+l 
2N+K 



). If / C 



0. 



K+l 

2N+K 



we are done. Otherwise, we need 

□ 



another shift and / (t) := f\ (t — a) will do, where / = (a, b) . 

Here is a simple but slightly disappointing corollary to theorem |21 
Corollary. M and D may be significantly different even when they are both small. 
Proof. Take b — K and A = K 2 . Then easily 

M(S) = + 0(R- 2 ) D(S) = (log 2) AT 1 + 0(A~ 3 / 2 ). □ 
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